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HIGHLIGHTS 


►  This  article  quantifies  and  models  the  impact  of  catalyst  layer  area  loss  during  operation  on  fuel-cell  performance. 

►  Oxygen  transport  through  the  ionomer  film  is  shown  to  become  limiting  due  to  the  reduced  area  for  reaction. 

►  Both  modeling  and  experiments  are  undergone. 
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During  proton-exchange-membrane-fuel-cell  degradation  various  changes  occur  in  the  membrane- 
electrode  assembly,  and  consequently  the  performance  decays.  Therefore,  the  correlation  between 
physical  property  changes  and  the  performance  decay  is  important.  In  this  paper,  we  focus  on  this 
correlation  and  have  developed  a  quantitative  model  to  estimate  the  performance  decay  after  degra¬ 
dation.  Comparison  of  electrochemical  surface  area  (ECSA)  and  polarization  curves  before  and  after 
potential  cycling  demonstrates  that  there  was  more  performance  decay  than  that  estimated  due  only  to 
ECSA  decrease,  especially  at  higher  current  densities.  This  large  performance  decay  is  estimated  and 
quantitatively  simulated  to  be  due  to  the  increase  in  oxygen  transport  resistance  caused  by  the  oxygen 
transport  through  the  ionomer  film  coating  the  Pt/C  particles  in  the  cathode  catalyst  layer.  It  is  shown 
that  the  data  can  be  fit  by  assuming  that  the  effective  ionomer  film  area  is  proportional  to  ECSA. 
According  to  our  results,  the  transport  resistance  increased  significantly  during  potential  cycling  and 
more  than  half  of  the  performance  decay  at  high  current  densities  originated  from  the  increased 
transport  overpotential.  The  effect  of  proton  overpotential  in  the  cathode  catalyst  layer  was  also  studied 
with  additional  performance  decay  witnessed  under  low-humidity  conditions. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Proton-exchange-membrane  fuel  cells  (PEMFCs)  are  one  of  the 
strongest  candidates  for  future  automobile  applications  due  to  their 
high  efficiency  and  clean  emissions.  In  automobile  applications, 
a  PEMFC  must  be  designed  not  only  to  work  under  a  wide  range  of 
conditions  but  also  to  maintain  its  performance  for  long-term  use 
over  years.  Therefore,  understanding  and  optimizing  the  degrada¬ 
tion  mechanisms  is  critical  throughout  the  entire  lifespan  of  the 
cell;  modeling  is  ideally  suited  to  examine  these  tradeoffs.  Here, 
‘degradation’  includes  both  physical  property  changes  and  perfor¬ 
mance  decay.  Typical  physical  property  changes  observed  in 
PEMFCs  include  decreasing  electrochemical  surface  area  (ECSA)  by 
dissolution  of  Pt,  carbon  corrosion,  thinning  of  the  catalyst  layer 
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and  membrane,  sintering  of  Pt  particles,  etc.,  while  performance 
decay  means  the  decline  of  the  polarization  curve  performance  due 
to  the  above  physical  property  changes.  There  have  been  plenty  of 
papers  published  recently  concerning  property  and  phenomeno¬ 
logical  changes  during  degradation  [1-13],  however,  the  correla¬ 
tion  between  these  physical  properties  changes  and  performance 
decay  has  not  been  fully  explained. 

Among  various  physical  property  changes,  ECSA  decrease  by  Pt 
dissolution  is  one  of  the  main  phenomena.  Assuming  only  a  decline 
in  activity,  the  cell  voltage  should  drop  by  the  same  degree  over  the 
entire  current-density  range  due  to  ECSA  decrease  based  on  the 
Tafel  equation 


RT . 

Vorr  -  — 


*ORR 


acF  Upti£ 


ref 


,ORR  J 


(1) 


where  tjorr  is  the  activation  overpotential,  R  is  the  ideal-gas 
constant,  T  is  the  absolute  temperature,  F  is  Faraday’s  constant, 
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Fig.  1.  Schematic  illustration  of  the  one-dimensional  PEMFC  model. 


^PtiOref.ORR  is  the  effective  exchange  current  density,  ac  is  the 
cathodic  transfer  coefficient,  and  i'orr  is  the  reaction  current 
density.  However,  significant  performance  decay  in  the  high- 
current-density  range  is  often  observed  in  experiments  that 
deviate  from  the  expected  decline  given  by  the  Tafel  equation 
[11-13].  Some  reports  attribute  this  to  ionomer  degradation  [14], 
but  a  quantitative  analysis  has  not  been  conducted. 

In  our  recent  study  [15,16],  the  oxygen  transport  resistance  of 
a  PEMFC  was  quantified  by  measuring  the  limiting  current  density 
using  diluted  oxygen  by  N2  or  He.  We  found  that  the  oxygen 
permeation  through  the  ionomer  film  is  a  major  portion  of  the  total 
transport  resistance  through  the  cathode  catalyst  layer,  especially 
for  low-Pt-loaded  catalyst  layers.  This  effect  has  also  been  shown  in 
various  modeling  studies  to  be  a  critical  issue  [17,18].  Therefore,  we 
anticipate  that  the  performance  decay  might  have  a  similar  effect 
with  the  oxygen  transport  resistance. 


The  aim  of  this  study  is  to  analyze  quantitatively  the  perfor¬ 
mance  decay  and  develop  an  estimation  technique  of  polarization 
curves  after  degradation,  especially  at  high  current  densities, 
because  such  a  technique  has  not  been  effectively  established  yet. 
To  analyze  the  mechanism,  we  use  both  experiments  and  a  ID 
mathematical  simulation  with  a  detailed  catalyst-layer  model  to 
evaluate  the  measure  polarization  curve.  Furthermore,  we  simulate 
the  impact  of  different  operating  conditions  on  the  performance 
after  degradation  such  as  under  low  humidity  and  low  oxygen  to 
ascertain  their  effects. 

2.  Modeling 

2.1.  One-dimensional  cell  model 

The  one-dimensional  model  of  each  component  of  the  PEMFC  ( see 
Fig.  1),  i.e.  gas-diffusion  layer  (GDL)  and  membrane-electrode¬ 
assembly  (MEA)  composed  of  catalyst  layers  (CLs)  and  membrane 
(Mem),  was  constructed  with  general  physical  equations  (mass 
balance,  heat  balance,  reactions,  etc.)  [19].  The  five-layer  ID  model 
includes  about  80  physical  properties  of  the  various  cell  components. 
The  governing  equations  of  each  PEMFC  component  are  summarized 
in  Table  1.  Additional  kinetic  and  physical-property  equations  are 
listed  in  Table  2.  Based  on  these  equations,  the  ID  model  of  each 
component  of  PEMFC  was  constructed  using  a  commercial  process 
modeling  software,  gPROMS®  [20].  The  components  were  connected 
to  each  other  by  defining  continuity  of  all  variables.  The  other  model 
parameters,  physical  properties,  are  in  Table  3.  These  models  are 
developed  with  the  following  assumptions: 

a.  Ideal  gas 

b.  No  liquid  water  in  the  cell;  excess  water  exists  as  supersatu¬ 
rated  vapor  or  mist. 


Table  1 

Governing  equations  for  steady  state  calculation. 


GDL 

CL 

Mem 

Energy  balance 

0  =  V-  (/<:qVT)  +  Q. 

Q  =  i2/Oe 

0  =  V-(fc,VT)  +  a 

(a"0da)Q  =  iH0R,H0R  +  (FiVH+)2/<  +  i>e  +  42000 fcM 

(cathode)  Q  =  i0RR(,0RR  +  0.25)  +  (F)VH02/<f  +  l2/^e  +  42000JGM 

0  =  V'(/<qVT)  +  Q 

Q  =  (fNH*)2/KM 

vG  =  -kG//iVP 

=  -kG//JLVP 

0  =  VNi 

(anode)  o  =  VNh2  +  i'hOr/2F 
(cathode)  Q  =  v.No;  +  i'orr/4F 

Ni  =  t^VPi 

Mass  balance  (gas  phase) 

VQ  E  (Cji  -  CJi/CtotalDff ) 

o  =  V'% 

0  —  V-(VVap  +7gm 

0  =  VNi 

(Crossover  gas) 

N[  =Ji+CiVG 

Ji  =  O^nd^Q 

Ni  =Ji+Qvc 

(ionomer  phase) 

- 

0  -  — V-Nwat; M  +Jgm 

^wat,M  =  ;  ^ll  —  ^wat  VCwa[M 

0  =  V-Nwat]yi 

^wat.M  =  ^H+  —  t^wat  VCwat[y] 

(Phase  change  flux) 

- 

Jgm  =  f<GM<Cvap  -  C/ap) 

- 

Charge  balance  (carbon  phase) 

0  =  V-(<TeV#c) 

(anode) Q  =  V  •  (<Te V#c)  -  iHOR 
(cathode)  0  =  V- (<TeV<l>c)  +  i'orr 

- 

(anode)  o  =  V  •  (FNH*  )  +  iH0R 

0  =  V-Nh+ 

(ionomer  phase) 

(cathode)  0  =  V(FNH*)-i0RR 

=  -CV#M 

FNHa=-KMV#M-* 
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Table  2 

Kinetic  equations  and  physical  property  calculation. 


HOR  kinetic  equations 

A  ■  (p h2  f 

*HOR  =  ^Pt*0  I  ^Fef exp  l^-^rr/HOR 

HOR  overpotential 

VHOR  =  -&C  +  &M-  ^HOR  [V] 

HOR  equilibrium  potential 


-  exp 


-acF_ 
RT  ~ 


-Vhor 


[Am 


Uhor  =  §l0§ 


[V] 


HOR  exchange  current  densities 
/ 2001 

V  R  l  Zef~T 


A  .  [20000  f  1  T 

AWo,hor  =  1-0  x  1011  exp  I  — „  I  = - = 


[A  m 


ORR  equilibrium  potential 
rrn  RT, 


Ut 


2  F 


;log(aw)  [V] 


ORR  standard  potential 
HqRR  =  1.23  -  0.0009(7  -  298.15)  [V] 

ORR  exchange  current  densities 


Ati$,RR  =  8.0  X  106 -exp ^38000  [A  m 

Oxygen  permeability  of  ionomer  :  measured  with  a  membrane  form 
W0l  =  5.0  x  10“15exp(1.4aw)-exp 


17000 

R 


[mol  s_1m_1Pa_1] 


Ionomer  water  content 
A  =  1.81  exp(1.89aw)  [-] 

Proton  conductivity  of  ionomer  :  measured  with  a  membrane  form 
12000  /  1  l\\  rc  _tl 

lSm  1 

Effective  proton  conductivity  of  ionomer  in  catalyst  layer 
[Sm-1] 


km  =  14.9aJ,75exp 


C„  - 

Kivr  —  ~ 
TM 


^eff 


Water  diffusivity  of  ionomer  :  measured  with  a  membrane  form 
Dw  =  (3.76  +  8.82aw  +  53.5 a2,  -  76.2a3,  +  28a4 ) 

/ 20000  /  1  I"* 


I ii  / 20000  /  1  l\\  r 

x10  exp(—  (^-rjj[n 


Effective  water  diffusivity  of  ionomer  in  catalyst  layer 

Dwat  =  ^MA/vat  [mV1] 

Electro  -  osmotic  drag  coefficient  of  ionomer 

£  =  1 

Binary  diffusivity  coefficient 


ff  _  c0  o.ooi. r1-75 

U  ^0P(,1/3+  l/3}2 


[m2s“ 


Knudsen  diffusivity  coefficient 


neff  _  fOf 
^Kundil 


/8RT  r  2  _ 1 1 

—rr  nrHs  1 

/  irMi 


c.  Convective  heat  transfer  is  ignored,  only  conduction  is 
considered. 

d.  Regarding  gas  diffusion,  Stefan-Maxwell  equations  in  GDL  and 
Knudsen  diffusion  in  CL  are  calculated,  and  other  diffusion 
effects  are  ignored. 

e.  Two  layered  GDL  is  treated  as  a  one-layer  GDL  of  average 
physical  properties. 

2.2.  Reaction  model 

An  agglomerate  approach  is  used  for  the  reaction  model  in 
the  catalyst  layer.  In  this  fashion,  the  oxygen-reduction-reaction 
(ORR)  kinetic  equations  with  ionomer  film  over  a  secondary 
carbon  (w/Pt)  agglomerate  for  PEMFCs  have  been  proposed 
recently  [17,18,21—29].  In  this  formulation,  the  transfer  current 
density  is  described  as 


ioRR  =4 

1  02 


+  : 


(2) 


{<^02, ion  WV 


where  P02  is  the  oxygen  partial  pressure  of  gas  phase  in  cathode  CL, 
Pq|  is  the  reference  oxygen  partial  pressure,  af*fn  is  the  effective 
ionomer  film  area,  <5ion  is  the  ionomer  film  thickness,  W02, ion  is  the 
oxygen  permeability  of  ionomer  and  /<orr  is  the  reaction  rate. 
Assuming  spherical  agglomerate  particles,  the  agglomerate  effect 
0ag g  is  expressed  as 


#agg  =  -2  (0COth0  -  1) 

<P 

where  0  is  Thiele  modulus, 


m  1 


KQRR 


.eff 

aSS\l  prefm  . 

1  02  v  02, ion 


(3) 


(4) 


and  rff  is  the  effective  agglomerate  radius.  The  kinetic  rate  equa¬ 


tion  for  the  ORR  is  given  by  a  Tafel  expression, 

/<ORR  =  0  -  ^PtOx)^PtioeORRexP  (  -  -^jVoRR 


(5) 


where  Apt  is  the  Pt  surface  area  per  unit  geometric  MEA  area.  The 
overpotential,  ijorr,  is  defined  as 


VORR  =  -^C  +  ^ion  +  ^ORR 


(6) 


Table  3 

Physical  properties  of  components  for  MEA1. 


Property 

Value  sources 

aGDL 

aCL 

Mem 

cCL 

cGDL 

unit 

Thickness 

measured 

5 

2  x  1(T4 

3  x  10"6 

3  x  lO"5 

5  x  10"6 

2  x  10"4 

m 

Porosity 

calculated 

£0 

0.75 

0.6 

- 

0.62 

0.75 

Tortuosity 

fitted 

TO 

3 

4 

- 

4 

3 

Electrical  conductivity 

measured 

CJe 

250 

60 

- 

60 

250 

S  ITT1 

Thermal  conductivity 

measured 

kq 

0.33 

0.5 

0.26 

0.5 

0.33 

W  m  1  K1 

Gas  permeability 

measured 

/<G 

1.2  x  10"13 

5  x  10-14 

- 

5  x  lO"14 

1.2  x  10"13 

m2 

Knudsen  radius 

assumed 

Tkund 

- 

4  x  lO"8 

- 

4  x  lO"8 

- 

m 

Ionomer  fraction 

calculated 

£m 

- 

0.2 

- 

0.15 

- 

Ionomer  tortuosity 

measured 

Tm 

- 

1 

- 

7 

- 

Transfer  coefficient 

assumed 

aa 

- 

1 

- 

- 

- 

Transfer  coefficient 

assumed 

0te 

- 

- 

- 

0.5 

- 

Phase  change  rate  coefficient 

assumed 

/<GM 

- 

4.0  x  105 

- 

4.0  x  105 

- 

s-1 

Effective  agglomerate  radius 

fitted 

reff 
'  agg 

- 

- 

- 

1.2  x  10-7 

- 

m 

Ionomer  film  thickness 

fitted 

5  ion 

- 

- 

- 

1  x  10~7 

- 

m 

Effective  ionomer  film  area 

fitted 

a?ff 

ion 

- 

- 

- 

1.5  x  105 

- 

m-1 
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where  Uorr  is  the  equilibrium  potential  and  0C  and  $jon  are  the 
potential  in  the  carbon  and  ionomer  phases,  respectively.  In  Eq.  (5), 
the  effect  of  platinum  oxides  (PtOx)  is  considered  in  the  0ptox  term, 
which  is  the  surface  coverage  ratio  by  PtOx  as  explained  below. 

In  cases  of  diffusion  limited  region  where  the  limiting  current 
density  (/nm)  is  observed  and  /<0rr  is  large  enough,  Eq.  (2)  can  be 
written  as  the  following  using  the  ionomer  related  transport 
resistance  (RCL,ion) 


*’orrI  k0] 


=  4FP( 


02 


R  =  4 F  (P02 \ 

CL>1°n  him  \RT  ) 


neff  m 
aion  ^  02, ion 


_  c  4  F 

—  l02,CL7 - 

him 


4FP, 


02 


1 


RTd  Rcui 


(7) 


(8) 


where  d  is  the  thickness  of  the  cathode  CL  and  C0 2, cl  is  the  oxygen 
concentration  of  gas  phase  in  the  cathode  CL. 

As  noted  in  Eq.  (5),  we  consider  the  effect  of  PtOx  coverage.  PtOx 
forms  on  the  surface  of  the  Pt  catalyst  above  0.6  V  in  the  presence  of 
water  or  oxygen.  These  oxides  disturb  and  inhibit  the  ORR  and 
subsequently  lower  the  cell  voltage.  Several  forms  of  oxide  and  its 
models  are  supposed  and  their  impact  have  been  measured 
[10,30,31].  Here,  a  simplified  model  to  simulate  PtOx  is  employed, 
which  assumes  a  single  electrochemical  reaction  governed  by 
a  Butler-Volmer  equation, 

d#ptox 


dt 


=  fc 


PtOx 


a'  F 

(1  -0ptox)exp(-j|jr?7ptox 


-  0ptoxexp 


-a'F 


RT 


-VptOx 


VPtOx  —  —&C  +  ^ion  +  hfptox 


(9) 


(10) 


where  /<ptox  is  the  reaction  rate  for  PtOx  formation,  a'a  and  a'c  is  the 
anodic  and  cathodic  transfer  coefficient  for  PtOx  formation, 
respectively.  Uptox  is  the  PtOx  equilibrium  potential.  From  Eq.  (9), 
the  PtOx  coverage  ratio  under  steady-state  conditions  can  be 
expressed  as 


exp 


^PtOx 


a'F 

RT 


VptOx 


exp 


RT 


VptOx 


exp 


-a'cF 

~irr 


VptOx 


The  calculated  relation  between  6ptox  and  Oc  -  (f>[on  is  shown  in 
Fig.  2.  Using  Eq.  (11)  in  Eq.  (5),  Eq.  (5)  can  be  rewritten  as  (see 
Appendix  A  for  the  detailed  derivation) 


^orr- 


RT 


-In 


ref  > 
ORR 


(a/a  +  a'c  +  ac)F  l  Iorr 
+  const  (Cell  voltage>>0. 86V)  (12) 

Therefore,  to  keep  the  Tafel  slope  nearly  70  mV/decade  at  353  K, 

<  +  <  +  aC  =  l  (13) 

Hence,  we  used  parameter  values  of  a'a  =  a'c  =  0.25,  ac  =  0.5, 

and  Uptox  =  0.76  V  from  literature  data  [30]. 

2.3.  Polarization-curve  fitting 

The  measured  polarization  curves  are  fitted  using  the  following 

methodology.  The  fitted  parameters  as  well  as  other  parameters  are 

shown  in  Table  3. 

(1)  Most  of  the  physical  properties  of  GDL  and  membrane  are 
assigned  from  existing  database.  For  example,  we  used 
measured  values  for  electrical  conductivity  and  thermal 
conductivity  of  GDL  with  direct  steady  state  method.  The 
contact  heat  resistance  was  included  in  the  bulk  value. 

(2)  Interfacial  electric  resistance  of  the  cell  is  fit  to  40  mQcnrr2 
from  the  measured  cell  resistance  with  the  current  interrupting 
method. 

(3)  Ionomer  tortuosity  is  calculated  from  the  proton  impedance 
measurement  of  CL. 

(4)  Parameters  of  the  CL  such  as  volume  fraction  of  pore  and  ion¬ 
omer  are  calculated  from  the  thickness  of  CL,  carbon/ionomer 
ratio,  and  Pt/carbon  ratio. 

(5)  Tortuosity  of  GDL  is  fit  to  the  transport  resistance  results. 
Tortuosity  of  CL  is  assumed  to  be  small  as  experimentally 
reported  [16]. 

(6)  Aptio, orr  x  d  is  fit  to  the  polarization  curves  (both  initial  and 
cycled)  at  the  low-current-density  region  (<0.2  A  cm-2). 

(7)  <5ion/af0fn  x  d  is  fit  to  the  polarization  curves  (at  initial  stage) 
with  1%  oxygen  concentration. 

(8)  Effective  agglomerate  radius  (rf^g)  is  obtained  by  fitting  the 
initial  polarization  curve  with  1%  dry-oxygen  concentration. 


(11)  3.  Experimental 


Electrodes  of  the  Pt/C  catalysts  were  fabricated  by  a  direct-spray 
method  onto  the  perfluorinated-polymer  type  membrane.  The 
anode  and  cathode  Pt  loadings  were  0.05  and  0.1  mg-Ptcm-2  (30  wt 
%  Pt/C  Tanaka  Kikinzoku  Kogyo),  respectively,  and  the  cathode 
ionomer/carbon  weight  ratio  (I/C  ratio)  was  0.75.  The  same  fabri¬ 
cation  method  was  used  for  other  MEAs  used  in  this  report  with 
different  species,  loadings,  and  Pt/C  weight  ratios. 

The  fabricated  MEA  (MEA1)  was  placed  in  a  13  cm2 
(3.6  x  3.6  cm)  test  cell  for  the  polarization,  cyclic-voltammetry 
(CV),  and  potential-cycling  measurements.  For  the  initial 
transport- resistance  measurement,  a  1  cm2  (lxl  cm)  test  cell  was 
used  to  ensure  high  stoichiometry  for  the  diluted  oxygen  feed. 
MEAs  were  sandwiched  by  GDLs  consisting  of  PTFE  treated  carbon 
paper  (200  pm  thickness)  and  a  microporous  layer  with  an 
assembly  pressure  of  1.0  MPa.  Porous-metal  gas  channels  (Ti  with 
gold  plating)  were  used  in  a  cross-flow  orientation. 

MEA1  was  tested  for  initial  H2/Air  polarization  performance  and 
CV  and  transport-resistance  measurements.  The  cell  conditions  are 
summarized  in  Table  4.  Polarization  curves  were  measured  under 
fully  humidified  conditions  with  a  cathodic  sweep  rate  of  50  mA  s-1 
using  a  commercial  fuel-cell  test  station  (850  series®,  Scribner 
Associates).  For  the  transport-resistance  measurements,  polarization 
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Table  4 

Cell  conditions. 


Atmosphere 
(flow  rate)  /  NL  min-1 

Anode  Cathode 

Pressure  / 
kPa  abs 

Temperature  / 
°C 

RH  / 

% 

CV  measurement 

H2  (2) 

Fluorinert™ 
FC-3283  (-) 

100 

40 

100 

Potential  cycling 

H2  (2) 

N2  (2) 

100 

80 

100 

Polarization  curve 

measurement 

H2  (2) 

Air  (2) 

200 

80 

90 

Transport  resistance 
measurement 

H2  (2) 

1%  02  99%  N2 
or  He  (2) 

150 

80 

90 

curves  were  measured  using  1  %  02/99%  N2  or  He  dry  gas  mole  fraction 
with  a  cathodic  sweep  rate  of  5  mV  s-1  using  an  AS-510-ECA® 
potentiostat  (NF  Corporation).  Unintentional  flooding  or  membrane 
dehydration  was  avoided  by  monitoring  the  high-frequency  cell- 
resistance  and  cell-voltage  stability. 

CV  measurements  were  performed  using  an  AS-510-ECA® 
potentiostat  (NF  Corporation)  by  sweeping  the  potential  between 
0  and  1.0  V  vs.  RHE  at  40,  50,  60  and  80  mV/s  sweep  rates.  While 
measuring  the  voltammogram,  humidified  H2  gas  was  fed  into  the 
anode  flow  field  and  the  cathode  electrode  was  filled  with  Fluo- 
rinert™  FC-3283  (perfluorohydrocarbon  liquid,  3M™)  [32].  If  N2  gas 
is  fed  into  cathode  flow  field,  then  hydrogen  evolution  would 
contribute  significantly  to  the  total  current  measured  at  the  low- 
potential  region  (<0.2  V)  of  the  hydrogen-adsorption  curves,  and 
this  effect  could  lead  to  the  overestimation  of  the  ECSA  [33].  By 
adopting  the  above  method  instead  of  using  N2  gas,  it  is  possible  to 
make  the  hydrogen-evolution  potential  lower  and  mitigate  such 
error.  The  ECSA  was  determined  from  the  total  charge  for  hydrogen 
adsorption  by  integrating  between  a  straight  baseline  drawn  from 
the  double-layer-capacitance  region  (0.4  V)  and  the  final  minimum 
potential  (0  V)  in  the  cathodic  CV  spectrum.  The  total  charge  for 
hydrogen  adsorption  depends  on  the  sweep  rate  and  the  effect  of 
hydrogen  evolution  on  the  total  charge  decreased  as  the  sweep  rate 
increased  (Fig.  3).  Therefore,  extrapolation  to  infinite  sweep  rate 
was  used  to  determine  the  charge;  this  method  can  minimize  the 
effect  of  hydrogen  evolution  [33].  Then,  210  pC  cm-2-Pt  capacitance 
was  used  to  determine  the  ECSA. 

Finally,  MEA-voltage-cycling  experiments  were  conducted  up  to 
49000  cycles  using  a  square-wave  potential  from  0.7  V  to  0.9  V  vs 
RHE  for  4  s  each  under  fully  humidified  conditions  with  H2/N2  (see 
Table  4  for  other  conditions).  During  the  cycles,  CV  spectra  and 
polarization  curves  were  measured  after  5000,  13000,  28000, 
39000  and  49000  cycles. 


Fig.  3.  Total  charge  for  hydrogen  adsorption  as  a  function  of  potential  sweep  rate  with 
different  sweep  rates  of  40,  50,  60  and  80  mV/s  for  MEA1  at  initial  stage. 


4.  Results  and  discussion 

4.1.  Transport  resistance  in  the  catalyst  layer 

The  initial  polarization  performance  is  shown  in  Fig.  4.  The  tests 
used  1%  02  in  N2  or  He  balance  gas  and  the  limiting  current  density 
was  0.22  Acm-2  and  0.32  Acm-2  for  the  N2  and  He  cases,  respec¬ 
tively.  The  total  transport  resistance  in  the  cathode  (Rtot)  was  then 
calculated  from  the  limiting  current  and  the  oxygen  concentration 
in  cathode  flow  field  using 

Ktot  =  Q)2,Gc4F//iim  =  RcL,tot  +  FGDL  (14) 

where  Co2,gc  is  the  oxygen  concentration  of  gas  phase  in  cathode  flow 
field.  To  distinguish  the  transport  resistance  in  CL  (Rc L,tot)  from  the 
transport  resistance  in  GDL  (Rgdl),  the  difference  of  molecular 
diffusivities  is  used  for  the  two  balance  gases  (D02  in  He/Do2  in  N2  =  2.19 
at  80  °C  and  150  kPa).  As  the  average  pore  size  of  GDL  is  10-100  pm, 
most  diffusion  through  the  GDL  occurs  by  molecular  diffusion.  On  the 
other  hand,  Rcptot  consists  of  Knudsen  diffusion  (pore  size  of 
0.01-0.1  pm)  and  permeation  resistances  through  the  ionomer  film, 
which  are  independent  of  the  balance  gas.  Using  the  above  relations, 
Rcptot  was  calculated  as  28  s  itT1  and  Rgdl  was  calculated  as  37  s  m-1 
and  17  s  itt1  for  the  N2  and  He  conditions,  respectively. 

From  the  initial  ECSA  measurement,  the  area  was  calculated  to 
be  71  m2  m-2  at  the  initial  stage.  Because  this  value  almost  agrees 
with  the  total  Pt  surface  area  (APt)  estimated  from  the  average  Pt 
radius  of  4.1  nm  obtained  from  XRD  measurement,  the  whole  Pt 
surface  on  MEA1  is  considered  to  be  electrochemically  active. 
Fig.  5(a)  shows  the  measured  ECSA  during  the  potential  cycling  for 
MEA1,  and  it  is  found  that  the  ECSA  decreased  to  27  m2  m-2  at  the 
end  of  the  potential  cycling  test.  In  terms  of  polarization,  Fig.  5(b) 
shows  the  measured  curves  at  each  stage  including  both  the  data 
and  the  model-fit  lines.  From  the  curves,  the  performance  decayed 
over  the  entire  region.  The  measured  potential  difference  at  1.7 
Acirr2  at  the  end  of  the  test  was  170  mV  which  is  more  than  5  times 
larger  than  that  at  0.2  Acm-2.  In  the  experiments,  the  large 
performance  decay  at  high  current  density  (>1  Acm-2)  was 
observed  as  well  as  the  small  performance  decay  at  kinetic  region. 
The  simulated  polarization  curves  described  the  same  tendency. 
The  performance  decay  both  at  0.2  and  1.7  Acm-2  are  shown  in 
Fig.  5(c).  The  performance  decline  slope  of  1.7  Acm-2  is  significantly 
steeper  than  that  of  0.2  Acm-2,  which  is  well  reproduced  by  the 
simulation. 

The  oxygen  transport  resistance  in  a  typical  PEMFC  consists  of 
molecular  diffusion  resistance  in  GDL,  Knudsen  diffusion  resistance 
in  the  gas  phase  of  the  catalyst  layer  (Ra,gas),  and  permeation 


Current  density  (A/cm2) 


Fig.  4.  Polarization  curves  of  limiting  current  measurement  of  MEA1  at  initial  stage  in 
1%  02  (dry)  with  balance  gases  of  N2  or  He,  at  80  °C,  150  kPa  and  90%  relative  humidity. 


S.  Jomori  et  al.  /  Journal  of  Power  Sources  215  (2012)  18-27 


23 


0  20000  40000  60000 

Cycles  (-) 


Current  density  (A/cm2) 


♦  0.2A/cnr  (EXP) 
a  1 ,7A/cnr  (EXP) 

—  0.2A/cnr  (SIM) 

—  1.7 A/cm2  (SIM) 


20000  40000 

Cycles  (-) 


Fig.  5.  (a)  Measured  ECSA  with  50  mV/s  (inset  of  (a))  CV  spectra,  (b)  Measured 
polarization  curves  of  MEA1  during  the  potential  cycling  and  model-fit  lines,  (c)  Cell 
voltage  with  respect  to  cycles  at  0.2  and  1.7  Acm-2. 


resistance  through  the  ionomer  film  surrounding  the  catalyst  layer 
agglomerates  (ftcuon)  [16].  The  total  transport  resistance  (Kcl, tot)  in 
the  CL  can  be  written  as 

^CL,tot  —  ^CL,ion  +  ^CL,gas  (15) 

The  schematic  illustration  of  the  effective  ionomer  film  area 
(aj^)  and  the  Pt  surface  area  (Apt)  is  shown  in  Fig.  6(a).  af^  is  the 
ionomer  film  area  surrounding  a  Pt  particle  that  is  effectively  used 
for  oxygen  permeation.  One  can  assume  that  af^n  is  roughly 
proportional  to  Apt 

afofn^Pt  (16) 

Accordingly,  using  Eq.  (16)  in  Eq.  (7),  Rc L,i0n  is  expressed  as 

^CL,ional/^Pt  (17) 

Combining  Eq.  (15)  with  Eq.  (17),  ftcutot  is  rewritten  as 


+  +  Q  (18) 

where  Cj  and  C2  are  constants.  This  relationship  indicates  Rc L,tot  is 
inverse  proportional  to  ECSA,  provided  the  ionomer  film  thickness 
is  fixed.  On  the  other  hand,  the  experimental  results  of  the  relation 


Fig.  6.  Schematic  illustration  of  Pt  surface  area  (Apt),  effective  ionomer  area  (af^),  and 
ionomer  film  thickness  (<5Ion)  in  catalyst  layer,  (a)  At  initial  stage,  (b)  After  degradation. 


between  ECSA  and  l?cL,tot  for  MEA1  are  shown  in  Fig.  7.  The  other 
experimental  results  of  various  MEAs  with  different  Pt  species, 
loadings  and  Pt/C  weight  ratios  are  also  shown  in  the  same  figure. 
The  dotted  line  is  the  fitting  curve  using  Eq.  (18)  (C?  =  1700  s  m_1 
and  C2  =  5  s  m-1),  which  shows  good  agreement,  thereby  indi¬ 
cating  the  assumption  of  Eq.  (16)  is  rational. 

ftcctot  can  be  calculated  using  Eq.  (18)  from  the  measured  ECSA 
during  the  potential  cycling.  The  calculated  Rc L,tot  and  ratio  of  af^fn 
change  are  shown  in  Fig.  8.  Rc L,tot  increased  from  28  s  m”1  at  initial 
stage  to  67  s  m_1  as  af^fn  decreased  by  62%.  The  explanation  of 
increased  ftcutot  after  potential  cycling  can  be  explained  as  in 
Fig.  6(b).  Pt  particles  in  the  cathode  catalyst  layer  dissolve  at  high 
potential  and  either  diffuse  into  the  membrane  or  merge  onto 
larger  Pt  particles  (Ostwald-ripening),  resulting  in  the  total  ECSA 
decrease.  Therefore,  the  mechanism  of  this  increase  of  ftcutot  is 
similar  to  that  of  the  case  of  low  Pt  loading  [15-17],  that  is  to  say, 
a  higher  flux  through  the  ionomer  film  is  needed  to  maintain  the 
same  total  current  density  after  ECSA  decreases.  Thus,  the  oxygen- 
concentration  slope  within  the  ionomer  film  becomes  steep, 
resulting  in  lower  oxygen  concentrations  at  the  surface  of  Pt 
particles  (see  Eq.  (2)). 
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Fig.  7.  The  relation  between  ECSA  and  the  transport  resistance  in  cathode  CL  at  the 
initial  stage  of  MEA1  ( ♦ ),  MEAs  with  different  Pt  species,  loadings  and  Pt/C  weight 
ratios  ( □ ),  and  the  fitted  line  of  Eq.  (18)  (-).  All  data  are  measured  at  80  °C,  150  kPa  and 
90%  relative  humidity. 
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Fig.  8.  The  calculated  ratio  of  change  of  af®,  (A)  and  the  transport  resistance  in 
cathode  CL  (RcL,tot)  (•)  of  MEA1  using  Eq.  (18)  at  80  °C,  150  kPa  and  90%  relative 
humidity. 


4.2.  Understanding  losses 

As  the  performance  decay  strongly  depends  on  current  density, 
the  breakdown  of  the  increased  overpotential  during  the  potential 
cycling  was  studied.  The  total  increased  overpotential  in  the 
cathode  CL  (=  total  performance  decay)  is  classified  into  three 
sources:  activity  and  oxygen  and  proton  transport.  The  activity 
overpotential  is  mainly  caused  by  the  ORR  activity  decrease; 
therefore  the  performance  decay  only  with  the  activity- 
overpotential  effect  is  calculated  by  reducing  APt  based  on  Eq.  (5); 
here  Fig.  5(a)  is  used  for  the  changing  ratio  of  APt.  The  increase  in 
the  oxygen  transport  overpotential  is  caused  by  the  increased 
oxygen-transport  resistance  in  the  cathode  CL  as  written  in  Eqs.  (7) 
and  (8).  Therefore,  the  increased  transport  overpotential  is 
expressed  by  the  gap  between  the  following  two  calculation 
results.  One  is  the  calculation  result  of  reducing  only  APt  and  the 
other  is  the  calculation  result  of  reducing  both  APt  and  a^n  using  the 
relation  of  Eq.  (18).  The  increase  in  the  proton-transport  over¬ 
potential  is  caused  by  changing  the  proton-transport  distance  in 
the  ionomer  by  shifting  the  reaction  distribution  across  the  cathode 
CL.  The  proton-transport  resistance  overpotential  (tjh.o.)  is  esti¬ 
mated  by  the  ionomer  potential  drop  from  the  membrane 
boundary  to  the  average  reaction  center  in  the  cathode  CL 

l 

J  IorrK^orr  (z)dz 

Vh.CL  =  ^ion  +  UoRR  ~  @C  ~  - - j -  (19) 

J  !orr(z)c1z 

0 

where  z  is  the  cross-sectional  MEA  dimension  from  the  membrane. 

The  polarization  curve  after  49000  cycles  was  analyzed  using 
the  above  method.  The  three  curves  (shown  in  Fig.  9(a))  were 
calculated:  initial  stage  [A],  only  increase  in  the  activity  over¬ 
potential  [B],  and  increases  in  both  the  activity  and  transport 
overpotentials  [C].  The  figure  indicates  that  most  of  the  perfor¬ 
mance  decay  less  than  1.0  AcirT2  is  caused  by  the  increase  in  the 
activity  overpotential,  while  over  1  Acm-2  the  transport  over¬ 
potential  increase  is  dominant  (see  Fig.  9(b)).  At  1.7  Acm-2  the 
transport  overpotential  increase  was  more  than  100  mV  compared 
to  5  mV  at  0.2  AcnrT2.  The  effect  of  the  increase  in  the  activity 
overpotential  was  about  40  mV  and  slightly  increased  at  1.7  Acm-2. 
The  increase  in  the  proton  overpotential  was  negligible  even  at  high 
current  densities  due  to  fully  humidified  inlets  (i.e.,  the  reaction 
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Fig.  9.  (a)  The  influence  of  each  overpotential  on  the  polarization  curve  after  49000 
cycles.  Dots  are  the  experimental  results  and  calculated  lines  are  overlaid.  Lines  are 
initial  [A],  only  with  activity  overpotential  increase  considered  [B],  and  with  transport 
overpotential  increase  added  [C].  (b)  The  breakdown  of  the  performance  decay  at  0.2 
and  1.7  Acm-2  after  49000  cycles. 
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distribution  within  the  cathode  CL  did  not  change  after 
degradation). 

4.2.1.  Limiting  case  1:  low-humidity  condition 

At  low-humidity  conditions,  however,  the  proton  overpotential 
can  become  significant.  To  study  this  effect,  simulations  were 
conducted  using  the  same  parameters  as  those  used  for  MEAl. 
Here  we  calculated  3  cases  under  low-humidity  conditions.  The 
first  case  is  the  initial  stage  [D],  the  second  is  after  degradation 
(49000  cycles)  with  uniform  ECSA  decrease  [E],  and  the  third  is 
after  degradation  with  nonuniform  ECSA  distribution  [F] 
(Fig.  10(a)  inset).  The  ECSA  distribution  across  the  cathode  CL 
becomes  nonuniform  after  potential  cycling  where  there  is  less 
ECSA  remaining  near  the  membrane  side  of  the  CL  [1,9,12]. 
Polarization  curves  were  then  calculated  for  each  case  (Fig.  10(a)). 
The  large  performance  decay  is  observed  in  case  [E]  and  [F],  and 
the  analysis  of  the  breakdown  of  the  overpotential  increase  was 
conducted  using  the  same  methodology  as  discussed  above. 
Fig.  11(b)  is  the  analysis  results  at  1  Acm-2  at  each  case  with 
different  humidity  conditions.  In  addition  to  the  transport  over¬ 
potential  increase,  there  are  some  additional  increases  in  the 
proton  overpotential. 

The  increased  proton  overpotential  is  explained  by  the  change  of 
the  reaction  distribution  across  the  CL  after  degradation.  In  general, 
under  low-humidity  conditions,  more  reaction  occurs  near  the 
membrane  side  of  the  CL  because  the  ionomer  proton-transport 
resistance  is  significant.  However,  in  case  [E],  the  reaction  distri¬ 
bution  is  unable  to  move  toward  the  membrane  side  because  of  the 
resultant  higher  oxygen-transport  resistance.  In  other  words,  the 
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Fig.  10.  (a)  Calculated  polarization  curves  of  MEA1  with  35%  relative  humidity  for  each 
case  of  [D_e],  [E_e]  and  [ F_e ] .  (inset  of  (a))  assumed  normalized  local  ECSA  across  the 
CL  for  MEA1  at  initial  stage  [D_e],  after  49000  cycles  with  uniform  ECSA  decrease  [E_e], 
and  after  49000  cycles  with  nonuniform  ECSA  decrease  [F_e].  (b)  Breakdown  of  the 
cause  of  the  performance  decay  at  1  Acm-2  at  each  case  for  different  humidity 
conditions,  (c)  Calculated  reaction  rate  distributions  across  cathode  CL  at  1  Acm-2  with 
35%  RH  for  each  case  of  [D_e],  [E_e]  and  [F_e]. 


Fig.  11.  (a)  Calculated  polarization  curves  at  initial  stage  and  after  degradation  with 
different  oxygen  concentration  under  90%  relative  humidity  at  80  °C  and  150  kPa. 
(b)  The  relation  between  increased  overpotential  at  1.0  Acm-2  and  dry  oxygen 
concentration  in  the  cathode  flow  field. 


the  average  oxygen  concentration  is  about  14%.  Generally,  the 
transport  overpotential  becomes  larger  as  the  current  density 
approaches  the  limiting  current  density  because  the  oxygen 
concentration  at  the  Pt  surface  approaches  zero.  It  is  important  to 
note  that  the  limiting  current  depends  on  humidity  as  well.  For 
example,  the  oxygen  partial  pressure  is  nearly  30%  larger  at  30% 
relative  humidity  compared  with  100%  relative  humidity  with 
150  kPa  total  pressure  at  80  °C.  Fig.  11(a)  shows  the  calculated 
polarization  curves  with  different  oxygen  concentrations  at  90% 
relative  humidity  at  80  °C.  Performance  decay  at  the  high-current- 
density  range  is  significant  with  low  oxygen  concentration. 
Fig.  11(b)  is  the  relation  between  the  average  dry  oxygen  concen¬ 
tration  in  the  cathode  flow  field  and  the  overpotential  increase 
after  degradation  at  1  cm-2.  It  is  evident  that  the  transport  over¬ 
potential  increase  becomes  significant  as  oxygen  concentration 
becomes  lower  while  the  activity  overpotential  increase  after 
degradation  is  almost  identical  and  independent  of  the  oxygen 
concentration. 


current  distribution  becomes  more  uniform  across  the  CL  (Fig.  10(c) 
[E_r])  in  order  to  minimize  the  total  potential  loss  comprised  of 
both  the  proton-  and  oxygen  transport  overpotentials.  This  current- 
distribution  change  is  amplified  in  case  [F],  where  there  is  less  ECSA 
remaining  near  the  membrane  side  (Fig.  10(c)[F_r]).  In  this  case,  the 
additional  proton  overpotential  occurs  because  the  proton- 
transport  distance  becomes  even  longer. 

4.2.2.  Limiting  case  2:  low-oxygen-concentration  condition 

In  practical  vehicle  use,  it  is  efficient  to  set  the  gas  stoichiom¬ 
etry  as  low  as  possible.  With  a  typical  air  stoichiometry  of  about 
1.4,  the  dry  oxygen  concentration  at  the  cell  outlet  becomes  7%,  so 


5.  Conclusions 

The  relation  between  physical  property  change  and  perfor¬ 
mance  decay  after  degradation  has  been  analyzed  using  a  detailed 
PEMFC  ID  model.  In  the  model,  the  effective  ionomer  film  area 
(afon),  which  describes  the  total  surface  area  of  ionomer  used  for 
effective  oxygen  permeation,  is  used  to  express  the  oxygen  trans¬ 
port  in  the  catalyst  layer  (CL).  The  analysis  of  polarization  curves 
after  potential  cycling  demonstrated  the  important  correlation  of 
afon  a  Apt  h  is  concluded  that  the  large  performance  decay  at  high 
current  density  is  mostly  explained  by  the  increased  transport 
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resistance  in  the  cathode  CL.  Therefore,  the  measurement  of  ECSA 
during  degradation  is  important  to  estimate  not  only  the  increase 
in  activity  overpotential  due  to  surface  area  loss,  but  also  the 
increase  in  the  transport  overpotential.  The  remaining  small  devi¬ 
ation  between  simulation  and  experiments  may  be  explained  by 
such  factors  as  catalyst-layer  structural  changes,  ionomer  degra¬ 
dation  (physical  property  changes),  and  carbon  corrosion.  For  the 
case  of  low-humidity  operation,  simulation  results  show  that 
the  performance  decay  expanded  after  degradation  compared  with 
the  high-humidity  condition.  This  effect  is  because  the  reaction 
current  distribution  shifts  away  from  the  membrane  interface  due 
to  the  oxygen-transport-resistance  increase  in  the  cathode  CL, 
which  is  compounded  by  the  nonuniform  ECSA  distribution  across 
the  CL.  By  taking  af^n  into  account  in  the  model,  the  performance 
decay  can  be  estimated  accurately. 
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Appendix  A 

Eq.  (12)  can  be  derived  as  follows.  Eq.  (5)  can  be  written  as 
RT 


VORR  ~ 


acF 


In 


^Ptl’aORR^ 


lORR 


+  ln(l  -  #ptox) 


(A-l) 


Substituting  Eq.  (11)  into  Eq.  (A-l),  tjorr  is  expressed  as 
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For  the  case  of  cell  voltage  higher  than  0.86  V,  where  the 
cathodic  reaction  is  negligible  compared  to  the  anodic  reaction, 
V ORR  is 
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By  using  the  relation  of  Eq.  (10)  into  Eq.  (A-3)  with  the  value  of 
Uorr  =  1-18  V  and  Uptox  =  0.76  V,  tiqrr  is 


RT  In  f Apt'<$>RR 1  (“a  +  <) 

VoRR 

Finally,  t^qrr  is  expressed  as 


(77orr  +  0.42)  (A-4) 


Nomenclature 

Ap  Pt  surface  area  per  unit  area  of  MEA,  m2  m-2 
aw  Water  activity,  - 

afon  Effective  surface  area  of  ionomer  film  on  catalyst  per  unit 
volume  of  CL,  m2  m~3 
C  Molar  concentration,  mol  m-3 

d  Catalyst  layer  thickness,  m 

Djj  Binary  diffusion  coefficient  of  i  and  j,  m2  s-1 
Dknud  Knudsen  diffusion  coefficient,  m2  s_1 
Dw  Diffusion  coefficient  of  water,  m2  s_1 
ECSA  Electrochemical  surface  area  per  unit  area  of  MEA,  m2  m-2 
F  Faraday  constant,  96485  C  mol-1 

io  Exchange  current  density,  A  m-3 

i  Current  density  per  unit  volume,  A  nrr3 

/iim  Limiting  current  density,  A  m  2 
J  Flow  rate,  mol  s-1 

Jgm  Evaporation/condensation  rate  between  gas  and  ionomer 

phase,  mol  m~3  s-1 

/<gm  Evaporation/condensation  coefficient  between  gas  and 

ionomer  phase,  s-1 
/<g  Gas  permeability,  m2 

/<orr  Reaction  rate,  mol  m-3  s-1 

kptox  Reaction  rate  for  platinum  oxidation,  mol  m~3  s_1 
kq  Thermal  conductivity,  W  m_1  K_1 

M  Molecular  weight,  — 

N  Flux,  mol  m-2  s-1 

Pi  Partial  pressure,  Pa 

P  Total  pressure,  Pa 

Q  Heat  generation,  J  m-3 

H<nud  Knudsen  diffusion  radius,  m 

rfffg  Effective  agglomerate  radius,  m 

R  Gas  constant,  J  moH1  K_1 

Ri  Transport  resistance,  s  m_1 

Rtot  Total  transport  resistance  in  CL,  s  m-1 
T  Absolute  temperature,  I< 

U  Equilibrium  potential,  V 

vG  Gas  velocity,  m  s_1 

z  Distribution  domain  for  cross-sectional  direction  of  MEA,  - 

Greek 

aa  Anodic  transfer  coefficient 

ac  Cathode  transfer  coefficient 

a'a  Anodic  transfer  coefficient  for  platinum  oxidation 

a'c  Cathode  transfer  coefficient  for  platinum  oxidation 

<5ion  Ionomer  film  thickness,  m 

7]  Activation  overpotential,  V 

e  Fraction,  - 

6  Pt  utilization  (agglomerate),  - 

#ptox  PtOx  coverage  ratio,  - 

k  Proton  conductivity,  S  m-1 

X  Water  content,  mol  mol-1 

Id  Viscosity,  Pa  s 

/dw  Chemical  potential  of  water  in  membrane,  J  mol-1 

?  Electro-osmotic  drag  coefficient,  - 

p  Bulk  density,  kg  mol-1  - 

<re  Electric  conductivity,  S  m-1 

t  Tortuosity,  - 

0  Theile  modulus,  - 

O  Potential,  V 

W02, ion  Oxygen  permeability  of  ionomer,  mol  s_1  m-2  Pa-1 
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Subscripts  and  superscripts 
CL  Catalyst  layer 

eff  Effective  value 
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gas  Gas  phase 

GC  Gas  channel 

GDL  Gas  diffusion  layer 

HOR  Hydrogen  oxidation  reaction 

i  generic  species 

ion  Ionomer  phase 

j  generic  species 

M  Membrane  phase 

ORR  Oxygen  reduction  reaction 

PtOx  Platinum  oxides 

ref  Reference  value 

vap  vapor  phase 

wat  water  phase 

0  gas  phase 
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